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Continued fraction expression of the Mathien series 

Xiaodong Cao*, Yoshio Tanigawa and Wenguang Zhai 


Abstract 

In this paper, we represent a continued fraction expression of Mathicu series by a contin¬ 
ued fraction formula of Ramanujan. As application, we obtain some new bounds for Mathieu 
series. 


1 Introduction 


The infinite series 

( 1 . 1 ) 


YA 2m 
2-^ <rn 2 + r 2 ) 2 ’ 

m —1 ' ' 


(r > 0 ), 


is called a Mathieu series. It was introduced and studied by Emile Leonard Mathieu in his 
book [T8] devoted to the elasticity of solid bodies. Since its introduction the series S(r) and 
its various generalizations have attracted many researchers, who established some remarkable 
properties of these series including the various integral representations, the asymptotic expan¬ 
sions, lower and upper estimates, see e.g. Cerone and Lenard [Sj, Frontczak m, Milovanovic 
and Pogany m, Pogany, Srivastava and Tomovski |21j, and references therein. 
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An integral representation for the Mathieu series CEO is given by Emersleben 


as 


( 1 . 2 ) 


1 f c 

5(r) = r / 

r Jo 


x 


■ sin (rx)dx. 


. P X 1 

/ o e 1 

The integral representation was used by Elbert HH to derive the asymptotic expansion of S(r ): 

OO t-, i 


(1.3) S(r)=J2(-Vr -^ = (r-Mx,), 

m=0 


where B 2n denote the even indexed Bernoulli numbers defined by the generating function 

y,n 


(1.4) 


X 


— 


t = X>A’ 

1 t—* n\ 


\x\ < 27r. 


n =0 


p/ /^\ 

Throughout the paper, we always use notation ip(z) = YTzf' Let (a n ) n >i and (6 n ) n >o be two 
sequences of real (or complex) numbers with a n 0 for all n € N. The generalized continued 
fraction 


(1.5) 


t = b 0 + 


a i 


b i + 


a 2 


= b 0 + 


CL 1 02 


h + b 2 + ■ 


= 6 0 + K 


=1 \ b 


^2+ * • 


( 1 . 6 ) 


= b 0 + K 


Ofc 


is defined as the limit of the nth approximant 

^4-n 

B n k =i 

as n tends to infinity. For the theory of continued fraction, see Cuyt, Petersen, Verdonk, 
Waadeland and Jones [9] or Lorentzen and Waadeland [16) . 

Let r > 0 and Rex > Let the continued fraction CF(r ; x) with a parameter r be defined 
by 

l 


(1.7) 

where for n > 1 

( 1 . 8 ) 


CF(r ; x) = 


(z-i) +A* 


Kn = - 


(® - \f + I (! + 4r2 ) + K “=i 

n 4 (n 2 + 4r 2 ) 1,9 9x 

4(2n — l)(2n + 1) ’ = ^(2n + 2n + 1 + 4r ). 


The main purpose of this paper is to establish the following continued fraction expression of the 
Mathieu series. 

Theorem 1. Let r > 0 and CL(r;x) be defined as (11.71) . For all positive integer k, we have 

k -1 


(1.9) 


S(r) = 27 


2m 


+ CF(v, k), 


—' (m 2 + r 2 ) 2 

m =1 v ' 

where the sum for k — 1 = 0 is stipulated to be zero. In particular, 
(1.10) S{r) = CF{r- 1). 
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2 Some preliminary lemmas 


In order to prove Theorem 1, we will prepare some lemmas. The following continued fraction 
formula of Ramanujan plays an important role in the proof of Theorem 1. 

Lemma 1 . Let x,l,m, and n denote complex numbers. Define 

(2.1) 

P = P(x, l, m, n ) 

T(|(x + / + m + n + l))r(^(a: + Z — m — ra + 1 )) T (^(x — l + m — n + 1)) T (|(x — l — m + n + 1)) 

T(!(x — l — m — n + 1)) T (^(x — l + m + n + l))T {^{x + l — m + n + 1)) T ^(x + l + m — n + 1)) 

Then if either l, m, or n is an integer or if Re x > 0, 


( 2 . 2 ) 


1 - P 
1 + P 


2 Imn 


;2 m 2 r 2 I 1 I K" 00 ( 4 (l * 2 -j 2 )(m 2 -j 2 )(n 2 -j 2 ) 

l -m -n + i + rv j=1 \( 2 j+i){x*-i*-rri*-ri*+ 2 p+ 2 j+i) 


Proof. This is Entry 35 of B. C. Berndt [5], p. 157, which was claimed first by Ramanujan [231 1 211 . 
The first published proof was provided by Watson [26]. For the full proof of Entry 35, we refer 
the reader to L. Lorentzen’s paper m- □ 

Lemma 2. The canonical contraction of bo + K (a n /b n ) with 


Ck = A 2k , D k = B 2k for k = 0 , 1 , 2 ,... 


exists if and only if b 2k 0 for k = 0 , 1 , 2 ,..., and is given by 

a\b 2 a 2 03 64 


OO 

bo + K 

n —1 \ Drj 


=b 0 +- 


040562^6 


T b\b 2 — 0364 + 62(04 + ^ 3 ^ 4 ) — 0566 + 64(06 + 6560 ) — ■ 

02n02n+l62n-262n+2 


— a, 2n+ ib 2n+2 + b 2n (a 2n+2 + b 2n+ ib 2n+2 ) - 


Proof. It follows from Theorem 12 and Eq. (2.4.3) of L. Lorentzen, H. Waadeland [H)], page 
83-84. For some applications, interested readers may refer to Berndt [5, p. 121, Eq. (14.2)] or 
[5, p. 157]. For a canonical contraction of a continued fraction and the related definitions, also 
see L. Lorentzen, H. Waadeland m, page 83. □ 

Lemma 3. 60 + K ( a n /b n ) « do + K {c n /d n ) if and only if there exists a sequence {r n } of 
complex numbers with ro = l,r n 7 ^ 0 for all n € N, such that 


(2.3) do = b 0 , c n = r n _ir n a n , d n = r n b n for all n € N. 

Proof. See Theorem 9 of L. Lorentzen, H. Waadeland [16], p. 73. □ 
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3 The proof of Theorem 1 

Lemma 4. Let r > 0 and Rex > then 


(3-1) 


CF{r ; x) — CF{r ; x + 1) = 


2 x 


(x 2 + r 2 ) 2 


Remark 1. In fact, Lemma 4 was guessed first by the multiple-correction method developed in 
mm- Mortici f20j made an important contribution in this direction. 


Proof. By applying Lemma 1 with ( x,l ) = (2x — l,2ri) and dividing both sides by 4 rmni, we 
obtain that for Re x > ^ 


(3.2) 

1 1 - P 

4ri 


imn(l + P ) i\2 , A r 2 m 2 „2 I 1 I yr°° ( 4(—4r 2 — j 2 ) (m 2 —j 2 )(n 2 —j 2 ) V 

v ' (Zx- i) +4r -m -n + 1 + JY j=1 ^ (2j+l)«2x-l)*+4r*-m*-n*+2p+2j+l)) 

Now let m tend to zero and n tend to zero, successively. On the right side, we arrive at 
(3.3) 1 


(0 T _ 02,1 I 4r 2 , K°° ( _ b'o 2 • 

ij + 1 + 4r + XV i j'=l \ i (2j+l)((2a;-l) 2 +2j 2 +2j+H-4r 2 ) 

On the other hand, from the definition of P, we see easily that lim m _>.oP = 1. A direct 
calculation with the use of L’Hospital’s rule gives 


(3.4) 


lim 


1 - P 


= lim 


1 


, 1 -P 1, 1 -P 1, d . 

„„ , ^ lim -= — lim - = — lim ——(1 — P) 

m ->0 m(l + P) m-t 0 1 + P m—>0 m 2 m—tO m 2 m -)0 dm 


1 ( . , n .. , ,n .. . . n .. , i 

=- j— 74 —- + x - n) + + x - m) + 74 -- + x + ri) — + x + m) j- 


By making use of L’Hospital’s rule again, and noting the following classical representation (e.g., 
see [1, Eq. 6.3.16, p. 259]) 


(3.5) 


^(* + l) = -7 + X) [u 


k =1 


1 


1 


k k + z 


{z + - 1,-2,-3, 


where 7 denotes Euler-Mascheroni constant, it follows from (13.41) that 


(3.6) 


lim lim 


1 -P 


= lim A f 


lim 


1 1 - P 


0 m—>o mn( 1 + P) n^o dn \m-> 0 ml + P 

=- (y/(x — ri) — ^(x + ri)) 


=5 E 


_ (x — ri + /c ) 2 ' (x + ri + fc ) 2 

fc—0 fc—0 / 


E 


= 2 ri 


x + k 


^ ((x + A ;) 2 + r 2 ) 


2P' 


fc =0 
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Combining (13.31) and (13.61) . we get that for Rex > ^ 


(3 7) V 2 ( x + k ) - 


jto (( x + k ) 2 + r2 ) 2 (2x - l) 2 + 1 + 4r 2 + j ( 


-4j 4 (j 2 +4r 2 ) 


(2j+l)((2z-l) 2 +2j 2 +2j+l+4r 2 ) 


=CF(r ; x). 


Here we used Lemma 3 in the last equality. It is not difficult to check that for Re x > \ 
(3.8) 



This completes the proof of Lemma 4. 

Proof of Theorem 1. We use the telescoping method. Theorem 1 follows readily from Lemma 
4. 

4 Some new inequalities for the Mathieu series 

The bounds for the Mathieu series attracted many mathematicians like Schroder |23], Emer- 
sleben [12] . Makai m and Diananda [10]. In the past twenty years, many authors like Alzer, 
Bagdasaryan, Brenner, Guo, Lampret, Milovanovic, Mortici, Pogany, Qi, Ruehr, Srivastava, To- 
movski, etc. have made important contributions to this research topic, see e.g Eua da na Gni 
EUE2] and references therein. Let us briefly recall some simple results. 

Mathieu [18J himself conjectured only the upper bound S(r ) < r -2 , r > 0, proved first by 
Berg [l]l. Makai m showed the double sided inequalities 


(4.1) 



Alzer et al. [2j improved the lower bound to 


(4.2) 



where the constant l/(2£(3)) and 1/6 are sharp. 

Milovanovic and Pogany m stated a composite upper bound of simple structure, 


(4.3) 



which is superior to (14.21) in the interval [0, y^(5 + 2\/3)/6 ~ 1.18772) 
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Let a\ = 1, b\ = (x — | + r 2 , for n > 1 


(4.4) 
and 

(4.5) 


«2n+l — 


n (n 2 + 4r 2 ) 
2 ( 2 n + 1 ) ’ 


n 


«2n — 


2 ( 2 n — 1 ) ’ 


b 2n +i = [x - 


1 


1 

- - + 


2 ) 4 2n + 1 ’ 

By Lemma 2, it is not difficult to prove that 


&2n — 1- 


(4.6) 


CF(r;x)=K(^). 


We let z = (x — tj ) 2 — 4 , ci = 2, di = 22 + 2r , 


(4.7) 


j c 2n +i = n(n 2 + 4 r 2 ), c 2n = n 3 , 

\ d 2n +i = 2 ( 2 n + 1)2 + 2 r 2 , d 2 n = 1 - 

It follows easily from Lemma 3 that 
Lemma 5. Let Rex > With the above notation, we have 


CF{r-x)= K 

71=1 \ On 


(4.8) 

Lemma 6. Assume x > For all positive integer l, then 

(4.9) 


K 1 / 

—i V d n 


2 1 


21-1 


K =- < CF(r-x) < K ^ 
n=1 \b n J n=1 \b n 

Proof. As the partial coefficients of the continued fraction K^Li j are positive, the assertion 
is deduced readily from the theory of the continued fraction and the first equality in Lemma 
5. □ 

The following theorem tells us how to obtain sharp bounds for the Mathieu series. 

Theorem 2. Let r > 0 and k,l € N. Let two sequences ( a(n)) n>1 , (b(n)) n>1 be defined by (|4.4I) 
and (14.5|) with x = k, respectively. Then 


(4.10) 


fc-i 

E 


2m 


( m 2 + r 2 ) 2 


21 /'n \ 

lkGtWxe 


2m 


21-1 


+ K — . 

(m 2 + r 2 ) 2 n=1 \b n 


In particular, 

(4.11) 

(4.12) 


+ 


1 


< S(r) < 


(l + r 2)2 5/2 + r 2 

2 , 4 
(1 + r 2 ) 2 + (4 + r 2 ) 2 + 13/2 + r 2 


+ 


1 


( 1 +r 2 ) 2 2 + r 2 

< 5(r) < 


+ 


+ 


(1 + ?’ 2 ) 2 (4 + r 2 ) 2 6 + r 2 
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Proof. (14.101) follows readily from Theorem 1 and Lemma 6. Taking (k,l) = (2,1) and (k,l) = 
(3,1) in (14.101) . respectively, we can obtain (14.111) and (14.121) . □ 

Remark 2. For comparison, our upper bound in (14.111) improves (4.3) when 0 < r < \J —2 + y/7 ~ 
0.803587. It is not hard to check that if 


r € 


-6 + 5C(3) 


'2 + x /-2 + llC(3) -5C 2 (3) 


2 + v /-2 + llC(3) -5C 2 (3)' 


C(3) - 1 


(0.0507096,4.44903), 


then our lower bound in (14.111) is superior to Alzer’s in (14.21) . In addition, the bounds in (14.121) 
are always superior to the bounds in (14.111) for all r > 0. 

Remark 3. Taking k = 2 in Theorem 1, letting r tend to zero, and then applying the second 
equality in Lemma 5, we can deduce the following continued fraction for Apery number £(3) 


,-von _ 1 i 1 l 3 l 3 2 3 2 3 

■ ) “ + 2' 2 ■ 1 + ~T~ + 2" • 3 + ~T~ + 2- • 5 +••• ' 


Also see Berndt [5], p. 155. 
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